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ABSTRACT 

We report on the results of a time-series photometric survey of NGC 2362, carried 
out using the CTIO 4m Blanco telescope and Mosaic-II detector as part of the Mon- 
itor project. Rotation periods were derived for 271 candidate cluster members over 
the mass range 0.1 < Af/M < 1.2. The rotation period distributions show a clear 
mass-dependent morphology, qualitatively similar to that in NGC 2264, as would be 
expected from the age of this cluster. Using models of angular momentum evolution, 
we show that angular momentum losses over the ~ 1 — 5 Myr age range appear to be 
needed in order to reproduce the evolution of the slowest rotators in the sample from 
the ONC to NGC 2362, as found by many previous studies. By incorporating Spitzer 
IRAC mid-IR measurements, we found that 3 — 4 objects showing mid-IR excesses 
indicative of the presence of circumstellar discs were all slow rotators, as would be 
expected in the disc regulation paradigm for early pre-main sequence angular mo- 
mentum evolution, but this result is not statistically significant at present, given the 
extremely limited sample size. 

Key words: open clusters and associations: individual: NGC 2362 - techniques: 
photometric - stars: pre-main-sequence - stars: rotation - surveys. 



1 INTRODUCTION 

NGC 2362 is a very young open clus ter, at a moderate dis - 
tance. We adopt the parameters of iMoitinho et all (|200ll ) 
for the remainder of this work: an age of 5~t 2 Myr, dis- 
tance 1480 pc and reddening E(B — V) = 0.10 mag. Note 
however that there is some controversy in the literature re- 
garding th e age of this cluster, with estimates ranging from 
3-9 Myr. iMavne et ail (|2007h have examined the relative 
ages of young open clusters, finding that NGC 2362 does 
indeed seem to be older than the ONC, but that it may be 
younger than IC 348, and they suggest an age of ~ 3 Myr. 

Nevertheless, the extremely young age, and lack of sub- 
stantial nebulosity or differential reddening in this cluster 
make it an ideal testing ground for th eories of pre-main se- 
quence stellar evolution. iDahml (|2005h carried out an exten- 
sive spectroscopic survey, concentrating on the Ha 6563 A 
and Li I 6707.8 A features, finding that a large fraction 
(~ 91 per cent) of the T Tauri star (TTS) population in 
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NGC 2362 is composed of weak-line emitters (WTTS), as 
expected given the canonical 5 Myr age for the cluster. They 
found that ~ 5 — 9 per cent of the TTS population are likely 
to still be undergoing accretion on the basis of the measured 
Ha equivalent widths, which is comparable to the fraction of 
objects found to e xhibit L-band exces s emission indicative 
of an inner disc bv lHaisch et alj (|200ll ). 



NGC 2362 has also been observed more recently 
with Spitzer using the IRAC (Infra-Red Array Cam- 
era) instrument in the 3.6 , 4.5, 5.8 and 8.0 fim bands. 
iDahm fc Hillenbrand! (|2007r t used these observations to de- 
rive spectral energy distributions for candidate NGC 2362 
members, and to characterise disc emission in objects ex- 
hibiting infra-red excesses. They found an upper limit to the 
fraction of objects exhibiting optically thick discs of ~ 7 ± 2 
per cent, with an additional ~ 12 ± 3 per cent of objects 
exhibiting weak or optically thin discs. 
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1.1 Evolution of stellar angular momentum 

A long-standing issue in our understanding of the evolution 
of stellar rotation rates on the pre-main sequence is that 
significant loss of angular momentum appears to be required 
in order to reproduce the observed rotation rates of low-mass 
stars. 

During the first few Myr of their evolution, low-mass 
pre-main sequence (PMS) stars contract by factors of ~ 2 — 
3. Hence, if angular momentum is conserved, all sub-solar 
mass stars, which are thought to have initial rotation periods 
of ~ 10 days, should rotate at periods shorter than 2 days 
after only 2 — 3 Myr. Observations at these ages show that 
some stars rotate much slower than this rate, and indeed 
the slowest rotators do not appear to spin up at all over the 
first few Myr of their evolution. 

The most popular scenarios for solving this problem in- 
voke additional angular momentum losses resulting fr om the 
inter action of the star with a circumstellar disc (e.g. iKonigll 
Il99ll) . In the "disc locking" scenario, this interaction re- 
sults in transfer of angular momentum from the star to the 
disc, which maintains a constant stellar rotation rate until 
the disc dissipates sufficiently to break the coupling. Alter- 
natively, the angular momentum may be lost through an 
enhanced stell ar wind, powered by th e accretion of material 
from the disc l|Matt fc Pudrit j|2005h . 

In both cases, the most basic prediction of the models 
is that there should be a correlation between rotation rate 
and the presence of discs, in the sense that the slow rotators 
should still be surrounded by discs, whereas for the rapid 
rotators, the discs should have dispersed. 

Much of the recent literature on rotation has been de- 
voted to the search for a correlation between rotation rate 
and disc indicators, including diagnostics of accretion (UV 
excess, spectral line diagnostics such as Ha, etc.) and of 
the thermal emission from the disc itself (infra-red excess). 
Many of the early studies were inconclusive, due to a com- 
bination of effects including small sample sizes, biases, de- 
pendence of the adopted disc indicators on other stellar pa- 
rameters (predominantly mass), and the use of ambiguous 
disc indicators. 

Two recent studies (|Rebull et al.ll2006l ; ICieza fc Balibed 

l2007h have claimed statistically significant detections of 
the expected correlation between slow rotat ors and cir- 
cums tellar discs, in the ONC (T ± 1 Myr; |Hillenbrandl 
Il997t ) and NGC 2264 (2 - 4 M yr; |Park et aIj|200Ch . respec- 
tively. Cieza & Baliber (2006) have also examined IC 348 
(2 — 3 Myr; iLuhman et al.l [20031 ). but the results were am- 
biguous due to small number statistics. These studies have 
all used the IRAC instrument on the Spitzer space telescope, 
which provides mid-IR photometry in four broad bands: 3.6, 
4.5, 5.8 and 8.0 (im. In particular, 8.0 infra-red flux ex- 
cesses provide a much more reliable indicator of the presence 
of circumstellar discs than shorter-wavelength IR excesses 
(e.g. in the K-band). There is thus mounting evidence in 
favour of the disc regulation paradigm, but nevertheless, it 
is important to verify that this result extends to older open 
clusters. 



1.2 The survey 

We have undertaken a photometric survey in NGC 2362 us- 
ing the CTIO 4m Blanco telescope and Mosaic-II imager. 
Our goals are two-fold: first, to study rotation periods for 
a sample of low-mass members, covering K and M spectral 
types, down to ~ 0.1 Mq, and second, to look for eclipsing 
binary systems containing low-mass stars, to obtain dynam- 
ical mass measurements in conjunction with radial velocities 
from follow-up spectroscopy. Such systems provide the most 
accurate determinations of fundamental stellar parameters 
(in particular, masses) for input to models of stellar evo- 
lution, which are poorly constrained in this age range. We 
defer discussion of our search for occultations to another 
publication (Miller et al., in preparation). 

These observations are part of a larger photometric 
monitoring survey of young open clusters over a range of 
ages and metalicities (the M onitor project: iHodgkin et al.l 
l2006l and lAigrain et alj|2007h . 

The remainder of the paper is structured as follows: the 
observations and data reduction are described in fJH and the 
colour magnitude diagram (CMD) of the cluster and candi- 
date membership selection are presented in Sj3] The method 
we use for obtaining ph otometric periods is summarised in 
(see| rwin et al 1 l2006l for a more detailed discussion). Our 
results are given in £j5j and f|6] summarises our conclusions. 

In several of the following sections, we use mass and ra- 
dius estimates for the cluster members. These were derived 
from the 5 Myr NextGen models, using the /-band magni- 
tudes, rather than V — I colour or V magnitude, for reasons 
discussed in H3.ll 



2 OBSERVATIONS AND DATA REDUCTION 

Photometric monitoring observations were obtained using 
the 4m CTIO Blanco telescope, with the Mosaic-II imager, 
on 8 x 1/2-nights between 2005 Dec 24 and 2006 Jan 06. 
This instrument provides a field of view of ~ 36' x 36' 
(0.37 sq.deg), using a mosaic of eight 2k x 4k pixel CCDs, 
at a scale of ~ 0.27"/pix. 

A single field in NGC 2362, centred on the star r CMa 
(see Figure [TJ, was observed for ~ 4 hours per night, in par- 
allel with another field in the cluster M50, the results from 
which will be published separately. Exposure times were 75 s 
in i-band, giving a cadence of ~ 6 minutes (composed of 2x 
75 s exposures plus 2x 100 s readout time, slewing between 
M50 and NGC 2362 during readout). We also obtained deep 
1^-band exposures (600 s, 450 s, 2 x 300 s and 150 s) which 
were stacked and used to produce a colour-magnitude dia- 
gram of the cluster. 

For a full desc ription of our data reduction steps, the 
reader is referred to llrwin et alj (|2007ah. Briefly, we used the 
pipeline for the INT wide-field survey l|lrwin fc LewislfeOQll ) 
for 2-D instrumental signature removal (bias correction, flat- 
fielding, defringing) and astrometric and photometric cali- 
bration. We then generated a master catalogue for each fil- 
ter by stacking 20 of the frames taken in the best conditions 
(seeing, sky brightness and transparency) and running the 
source detection software on the stacked image. The result- 
ing source positions were used to perform aperture photom- 
etry on all of the time-series images. We achieved a per data 
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Figure 1. Digitised sky survey (DSS) image of NGC 2362 cover- 
ing ~ 1.2° X 1° , showing the coverage of the present survey (8-chip 
mosai c tile), and the Spitzer IRAC field of lDahm fe Hillenbrand! 
1120071 1, showing the total coverage (larger box), and the smaller 
central region covered by all four bands. Crosses show the posi- 
tions o f the candidate cluster members from lDahm fc Hillenbrand! 
11200711 . 



also generated by the pipeline software l|lrwin fc LewislfeoOll ) 
to allow us to identify non-stellar or blended objects in the 
time-series photometry. 

Photometric calibratio n of our data w as carried out us- 
ing regular observations of iLandoltl (|l992l ) equatorial stan- 
dard star fields in the usual way. 

Light curves were extracted from the data for ~ 85 000 
objects, 56 000 of which had stellar morphological classifica- 
tions (~ 40 per cent of these are flagged as blended), using 
our standar d apert ure photometry techniques, described in 
I Irwin et al.l |2007al ). We fit a 2-D quadratic polynomial to 
the residuals in each frame (measured for each object as the 
difference between its magnitude on the frame in question 
and the median calculated across all frames) as a function 
of position, for each of the 8 CCDs separately. Subsequent 
removal of this function accounts for effects such as varying 
differential atmospheric extinction across each frame. Over 
a single CCD, the spatially- varying part of the correction 
remains small, typically ~ 0.02 mag peak-to-peak. The rea- 
sons for usi n g this technique are discussed in more detail in 
llrwin et al.l (1200731 ). 

For the production of deep CMDs, we stacked 20 i-band 
observations, taken in good seeing and photometric condi- 
tions, and all of the l/-band observations. The limiting mag- 
nitudes on these stacked images, measured as the approx- 
imate magnitude at which our catalogues are 50 per cent 
complete, were V ~ 24.4 and i ~ 23.6. 




i magnitude 



Figure 2. Plot of RMS scatter per data point (measured over 
the entire data-set of 8x 1/2- nights) as a function of magnitude 
for the i-band observations of a single field in NGC 2362, for 
all unblended objects with stellar morphological classifications. 
The diagonal dashed line shows the expected RMS from Poisson 
noise in the object, the diagonal dot-dashed line shows the RMS 
from sky noise in the photometric aperture, and the dotted line 
shows an additional 1.5 mmag contribution added in quadrature 
to account for systematic effects. The solid line shows the overall 
predicted RMS, combining these contributions. 

point photometric precision of ~ 2 — 4 mmag for the bright- 
est objects, with RMS scatter < 1 per cent for i < 19 (see 
Figure[2} . A signal-to- noise ratio of 5 (corresponding approx- 
imately to the detection limit for point sources on a single 
frame of the differential photometry) is reached at i ~ 22.7. 

Our source detection software flags any objects detected 
as having overlapping isophotes. This information is used, 
in conjunction with a morphological image classification flag 



3 SELECTION OF CANDIDATE LOW-MASS 
MEMBERS 

Catalogues of ca ndidate NGC 2362 members are available 
in the literature (|Moitinho et al.ll200ll ; lDahmll2005l ). but we 
elected to perform a new photometric selection using V ver- 
sus V — I CMDs from our data in order to properly match 
the relatively large field of view of our time-series observa- 
tions compared to these previous studies. 

3.1 The V versus V - I CMD 

Our CMD of NGC 2362 is shown in Figure O The V and 
i measurements were converted to the standard Johnson- 
Cousins photometric system using colour equations derived 
from our standard star observations: 

(V - I) = (V ccd - w) / 0.899 (1) 
V = V ccd + 0.005 (V - I) (2) 
/ = i cci ~ 0.096 {V - I) (3) 

Candidate cluster members were selected by defining 
an empirical cluster sequence 'by eye' to follow the clearly- 
visible cluster single-star sequence. The cuts were defined by 
moving this line along a vector perpendicular to the clus- 
ter sequence, by amounts k — a(V — I) and k + a(V — I) 
as measured along this vector, where a(V — /) is the pho- 
tometric error in the V — I colour. The values of k used 
were —0.25 mag for the lower line and 0.7 mag for the up- 
per line on the diagram, making the brighter region wider 
to avoid rejecting binary and multiple systems, which are 
overluminous for their colour compared to single stars. 1800 
candidate photometric members were selected, over the full 
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Figure 3. V versus V - I CMD of NGC 2362 from stacked im- 
ages, for all objects with stellar morphological classification. The 
cluster sequence is clearly visible on the right-hand side of the 
diagram. The boundaries of the region used to select photomet- 
ric candidate members are shown by the dashed lines (all objects 
between the dashed lines were selected). The reddening vector 
for Ay = 1.0 is shown at the right-hand side of the diagram. The 
mass scale is from the 5 Myr NextGen models llBaraffe et al ■Il998f) 
for M > 0.1 M Q , and the 5 Myr DUSTY models l|Chabrier et aD 
2000) for M < 0.1 Mq, using our empirical isochrone to convert 
the V magnitudes to I magnitudes, and subsequently obtaining 
the masses from these, due to known problems with the V magni- 
tudes from the models (see iJXTJl. The error bars at the left-hand 
side of the plot indicate the typical photometric error for an ob- 
ject on the cluster sequence. 



V magnitude range from V = 16 to 26, but the well-defined 
cluster sequence appears to terminate at M ~ 0.1 Mq, or 

V ~ 23.5, with a few candidate brown dwarfs found below 
this limit, but with high field contamination. 

We also consi dere d using the model i sochrones of 
iBaraffe et all (| 19981 ) and lChabrier et al l (|200Cl ) for selecting 
candidate members. The NextGen model isochrones were 
found to be unsuitable due to the known discrepancy be- 
tween these models and observations in the V — I colour 
for T cff < 3700 K (corresponding here to V - I > 2) This 
was examined in more detail by iBaraffe et alj |l998), and 
is due to a missing source of opacity at these temperatures, 
leading to overestimation of the V-band flux. Consequently, 
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Figure 4. Contamination estimated from Galactic models, mea- 
sured as the ratio of the calculated number of objects in each 
magnitude bin from the models, to the number of objects de- 
tected and classified as candidate cluster members in that magni- 
tude bin. Note that bins with contamination estimates > 100 per 
cent (where there were more objects in that bin from the Galactic 
model than were actually observed) have been truncated to 100 
per cent. 



when we have used the NextGen isochrones to determine 
model masses and radii for our objects, the /-band absolute 
magnitudes were used to perform the relevant look-up, since 
these are less susceptible to the missing source of opacity, 
and hence give more robust estimates. 



3.2 Contamination 

Since the cluster sequence is not well-separated from the 
field in the CMD, it is important to estimate the level of 
field star contamination in the sample of candidate cluster 
members. This has been done using two independent meth- 
ods. 



3.2.1 Galactic models 



The Besancon Galactic models ()Robin et all2003h were used 
to to generate a simulated catalogue of objects passing 
our selection criteria at the Galactic coordinates of NGC 
2362 (I = 238.2°, b = -5.5°), covering the total FoV of 
~ 0.35 sq.deg (including gaps between detectors). The same 
selection process as above for the cluster members was ap- 
plied to this catalogue to find the contaminant objects. A 
total of 1500 simulated objects passed these membership se- 
lection criteria, giving an overall contamination level of ~ 65 
per cent after correcting for bins where the number of ob- 
jects predicted by the models exceeded the number actually 
observed by ~ 20 per cent (we simply assumed 100 per cent 
field contamination in these bins). Figure [4] shows the con- 
tamination as a function of V magnitude. Note that this 
figure is somewhat uncertain due to the need to use Galac- 
tic models, and especially given the overestimation of the 
numbers of observed objects by the models, which lends low 
confidence to the results derived in this section. We there- 
fore pursue an alternative method for obtaining a limit to 
the contamination level below. 
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Figure 5. Upper limit to the field contamination level estimated 
from the radial distribution method (see text), plotted as a func- 
tion of V-band magnitude. 



3.2.2 Radial distribution 

An alternative method for generating a simple upper limit 
to the contamination level has also been applied to verify the 
results obtained from the Galactic models. Since the cluster 
is relatively centrally concentrated, we divided the popula- 
tion of candidate cluster members into two bins in radius, 
using as the centre the position of the star r CMa. Num- 
ber counts of objects selected using the CMD as candidate 
members were then produced inside and outside of a 15' ra- 
dius, chosen since it roughly divides the area covered into 
two, providing good statistics in both bins. By assuming the 
stars inside 15' are cluster members, and those outside are 
not, we can obtain a simple upper limit to the contamina- 
tion level by dividing these counts, and correcting for the 
relative area covered by each bin. 

Figure [5] shows the resulting distribution of contamina- 
tion as a function of magnitude, for comparison with the 
Galactic model results in Figure U The overall contami- 
nation level integrated over all magnitude bins from this 
method is 68 per cent, comparable to the estimate using the 
Galactic models from Section [3.2. II 



4 PERIOD DETECTION 
4.1 Method 

The method we use f or detection of pe riodic variables is 
described in detail in llrwin et al.l l|2006h . and we provide 
only a brief summary here. The method uses least-squares 
fitting of sine curves to the time series m(t) (in magnitudes) 
for all candidate cluster members, using the form: 



m(t) — rridc + ct sva(ust + ■ 



(4) 



where (the DC light curve level) , a (the amplitude) and 
cf> (the phase) are free parameters at each value of u over an 
equally-spaced grid of frequencies, corresponding to periods 
from 0.005 — 50 days for the present data-set. 

Periodic variable light curves were selected by evaluat- 
ing the change in reduced x 2 '- 

^Xf X(-* Xf, smooth ^ 0.4 

where xt is the reduced x 2 
respect to a constant model, and xt.s 




(5) 

of the original light curve with 
, is the reduced x 2 



hput period (days) 

Figure 7. Detected period as a function of actual (input) pe- 
riod for our simulations. Objects plotted with crosses had frac- 
tional period error < 10 per cent, open circles > 10 per cent. 
The straight lines represent equal input and output periods, and 
factors of 2, 3, 1/2 and 1/3. The curved lines are the loci of the 
±1 day" 1 aliases resulting from gaps during the day. The major- 
ity of the points fall on (or close to) the line of equal periods. 



of the light curve with the smoothed, phase-folded version 
subtracted. This threshold was used for the M34 data and 
appears to work well here too, carefully checked by exam- 
ining all the light curves for two of the detectors, chosen 
randomly. A total of 1008 objects were selected by this au- 
tomated part of the procedure. 

The selected light curves were examined by eye, to de- 
fine the final sample of periodic variables. A total of 271 
light curves were selected, with the remainder appearing 
non-variable or too ambiguous to be included. 



4.2 Simulations 

Monte Carlo simu lations were perf ormed following the 
method detailed in llrwin et a l. (2006), injecting simulated 
signals of 2 per cent amplitude and periods chosen following 
a uniform distribution on log 10 period from 0.1 to 20 days, 
into light curves covering a uniform distribution in mass, 
from 1.2 to 0.1 Mg. A total of 1015 objects were simulated. 

The results of the simulations are shown in Figure [6] as 
diagrams of completeness, reliability and contamination as 
a function of period and stellar mass. Broadly, our period 
detections are close to 100 per cent complete from 1.2 Mq 
down to 0.2 Mq, with remarkably little period dependence. 
Figure [7] shows a comparison of the detected periods with 
real periods for our simulated objects, indicating good reli- 
ability. 



4.3 Detection rate and reliability 

The locations of our detected periodic variable candidate 
cluster members on a V versus V — I CMD of NGC 2362 are 
shown in Figure [8] The diagram indicates that the majority 
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Figure 6. Results of the simulations for 0.02 mag amplitude expressed as percentages, plotted as a function of mass (left) and period 
(right). The simulated region covered 0.1 < M/Mq < 1.2 in order to be consistent with the NGC 2362 sample. Top panels: completeness 
as a function of real (input) period. Centre panels: Reliability of period determination, plotted as the fraction of objects with a given true 
period, detected with the correct period (defined as differing by < 20 per cent from the true period). Bottom panels: Contamination, 
plotted as the fraction of objects with a given detected period, having a true period differing by > 20 per cent from the detected value. 



of the detections lie on the single-star cluster sequence, as 
would be expected for rotation in cluster stars as opposed 
to, say, eclipsing binaries. 

The properties of all our rotation candidates are listed 
in Table [T] 



5 RESULTS 

5.1 NGC 2362 rotation periods 

Plots of period as a function of V — I colour and mass for the 
objects photometrically selected as possible cluster members 
are shown in Figure [9] Below ~ 0.7 M© (or M0), these dia- 
grams reveal a correlation between stellar mass (or spectral 
type) and the longest rotation period seen at that mass, with 
a clear lack of slow rotators at very low masses. This trend 
is also followed by the majority of the rotators in this mass 
range, with only a tail of faster rotators to ~ 0.5 day peri- 
ods, and very few objects rotating faster than this. This is 
very similar to wha t we found in the earlier NGC 2516 a nd 
NGC 2547 studies (llrwin et al.ll2007b| ; urwin et al.ll2007cl ). 

Above ~ 0.7 Mq, the distribution appears to show little 
mass dependence, with the slowest rotators all exhibiting 
~ 10 day periods. 

These morphological features do not appear to be a re- 
sult of sample biases. In particular, the simulations show 



that the survey is sensitive to much shorter periods than 
0.5 day, and the upper limit in detectable periods is not 
mass-dependent, so this cannot explain the sloping morphol- 
ogy of the upper envelope of rotation periods in Figure [9] at 
low masses. 

Figure [TU] may indicate a slight bias toward shorter pe- 
riods at low masses, if the distribution of amplitudes in the 
0.4 < M/Mq < 1.2 (top) and M < 0.4 M bins was the 
same, with no counterpart to the small group of objects at 
periods of ~ 5 — 10 days and amplitudes < 1 per cent visible 
in the lower panel. However, the result is not highly signif- 
icant due to small number statistics in the lower mass bin, 
and it is possible that these objects may be false positive de- 
tections, since at these amplitudes and periods, the effect of 
light curve systematics is sometimes difficult to distinguish 
from real variability. 



5.1.1 Period histograms 

In order to quantify the morphology of Figure [9] we have 
used histograms of the rotation period distributions in two 
broad mass bins, 0.4 < M/M < 1.2 and M < 0.4 M , 
shown in Figure [TT1 We have attempted to correct the distri- 
butions for the effects of incompleteness and (un)reliability 
using th e simulations descri bed in £1-4.21 following the method 
used in llrwin et al.l (|2006l ). The results of doing this are 
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Table 1. Properties of our 271 rotation candidates, including V and /-band magnitudes from our CCD photometry, Spitzer IRAC 
magnitudes in the 3.6, 4.5, 5.8 and 8.0 fim bands, where available, the peri od P, t-band ampli tude a; (in units of magnitudes, in the 
instrumental bandpass), interpolated mass and radius (from the models of iBaraffe et al. I ll998l, derived using the I magnitudes). Our 
identifiers are formed using a simple scheme of the cluster name, CCD number and a running count of stars in each CCD, concatenated 
with dashes. The full table is available in the electronic edition. Machine readable copies of the data tables from all the Monitor rotation 
period publications are also available at http://www.ast.cam.ac.uk/research/monitor/rotation/ 
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Figure 11. Period distributions for objects classified as possible photometric members, in two mass bins: 0.4 < M/Wq < 1.2 (upper 
row, corresponding roughly to K and early-M spectral types) and M < 0.4 Mq (lower row, late-M). The left-hand panels show the 
distributions plotted in linear period, and the right-hand panels show the same distributions plotted in log 10 period. The dashed lines 
show the measured period distributions, and the solid lines show the results of attempting to correct for incompleteness and reliability, 
as described in the text. 



shown in the solid histograms in Figure 1111 and the raw 
period distributions in the dashed histograms. 

The period distributions in the two mass bins of Fig- 
ure Qj] show clear differences, with the low-mass stars (M < 
0.4 Mq) showing a strongly peaked rotational period distri- 
bution, with a maximum at ~ 0.6 — 0.7 days, whereas the 
higher-mass stars (0.4 < M/Mg < 1.2) show a broader dis- 
tribution. We applied a two-sided Kolmogorov-Smirnov test 



to the corrected distributions to confirm the statistical sig- 
nificance of this result, finding a probability of 4x 10 -10 that 
the distributions were drawn from the same parent popula- 
tion. 

The implication of this result is that the observed mor- 
phology in Figure [9] and in particular the increase of the 
longest observed rotation period as a function of increasing 
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Figure 8. Magnified V versus V -I CMD of NGC 2362, for ob- 
jects with stellar morphological classification, as Figure [3] show- 
ing all 271 candidate cluster members with detected periods 
(black points). The dashed lines show the cuts used to select can- 
didate cluster members (see SI3.1I I. 



KO MO M3 MS 




Figure 9. Plots of rotation period as a function of dereddened 
V — I colour (top), and mass (bottom) for NGC 2362, deriving 
the masses us i ng th e 5 Myr NextGen mass-magnitude relations of 
iBaraffe et al,l lll998l ) and our measured /-band magnitudes. In the 
lower diagram, the greyscales show the completeness for 0.02 mag 
periodic variations from the simulations. 

clusters are close in "rotational age". Nevertheless, the ro- 
tation period results seem to indicate that the canonical age 
sequence, from youngest to oldest, of the ONC, NGC 2264 
and then NGC 2362, is correct. 



mass, a trend followed also by the bulk of the rotators, is 
real and statistically significant. 



5.2 Comparison with other data-sets 

5.2.1 Period versus mass diagrams 



Figure [12] shows a diagram of rotation period as a function 
of stellar mass for the ONC, NGC 2264 , NGC 2362, NGC 
2547 (~ 40 Mvr; iNavlor fc Jeffries! 120061) and NGC 2516 (~ 
150 Mvr; IJeffries. Thurston fc Hamblvll200ll ). Data sources 
for each cluster are indicated in the figure caption. 

The diagram clearly shows a gradual evolutionary se- 
quence, from a relatively flat mass- dependence of the rota- 
tion periods in the ONC (~ 1 Myr), to a sloping relation in 
NGC 2264 and NGC 2362, and an increasingly more pro- 
nounced slope at the older ages of NGC 2547 and NGC 
2516. 

The distribution for NGC 2362 appears to show a 
slightly more pronounced sloping relation than NGC 2264, 
but the distributions are very similar, indicating that these 



5.2.2 Comparison with rotational evolution models 

By comparing the observed rotation periods as a function of 
age with models of rotational evolution on the pre-main se- 
quence, we can constrain the amount of angular momentum 
loss required during this part of the evolution. 

In order to do t his, we have u sed a s imple model of solid 
body rotation from llrwin et al. I (|2007bl ), including angular 
momentum losses due to main sequence magnetically driven 
stellar winds (the effect of which is close to negligible during 
this age range), but specifically, not including any additional 
sources of angular momentum loss ("disc locking" or simi- 
lar) in order to see if the latter are required to explain the 
observations. Our models from the earlier NGC 2516 and 
NGC 2547 work incorporated a simple prescription for core- 
envelope decoupling, but at these early ages the stars are 
fully-convective, so it is not necessary to incorporate this 
effect here. 

Figure [TH] shows the results of attempting to evolve the 
measured rotation periods in the ONC forward in time to 
the age of NGC 2362 using these models. 
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Figure 10. Plot of amplitude as a function of period for NGC 
2362 in two mass bins: 0.4 < M/M < 1.2 (top) and M < 0.4 Mq 
(bottom). 



The model clearly cannot successfully reproduce the 
NGC 2362 distribution. Although there is some hint of the 
sloping relation at low masses visible in the evolved ONC 
distribution, the slowest rotators at M>0.6 M© are clearly 
rotating too rapidly. However, the lower envelope of fast ro- 
tators is reasonably well-reproduced throughout the mass 
range. 

We suspect that the mass determinations (which are 
based on /-band luminosity) are somewhat unreliable in the 
ONC, due to the presence of strong spatially- variable extinc- 
tion in this cluster. This introduces horizontal scatter in the 
period versus mass plane. It has been necessary to assume 
a typical (low) value of this parameter for many of the low- 
mass detections du e to a lack of suitab le estimates based on 
spectroscopy from iHillenbrandl (|l997r i. Hence, much of the 
scatter at low masses in the figure may be due to underesti- 
mation of the extinction, and hence of the masses, for highly 
reddened objects. 

However, at the high-mass end, the diagram is more re- 
liable, since the majority of the scatter due to extinction will 
only move objects to the right. Therefore, the result that the 
slowest rotators at these masses are rotating too rapidly is 
still significant. This indicates that an additional source of 
angular momentum loss is required, in order to hold the ro- 
tation periods of these stars approximately constant (since 
the slowest rotators in the ONC and NGC 2362 appear to 
have very similar periods of ~ 10 days). This is precisely 
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Figure 12. Rotation period as a function of stellar mass for 
(top to bottom): ONC, NGC 2264, NGC 2362, NGC 2547 and 
NGC 2516. Lower and upper limits (from usini data) are marked 
with arrows. The m asses were taken f rom the NextGen mass- 
magnitude relations l |Baraffe~et al. 1998) at the appropriate ages. 
The ONC data are fromlHerbst et al.1 j2002h ■ For NGC 2264 w e 
used the data of lLamm et al] lj2005T l and iMakidon et al.l l|2004h . 
The NGC 2547 a nd NGC 2516 rota t ion period data ar e from the 
Monitor project jlrwin et al]|2007bl : llrwin et al]|2007ch . 



the result found by many other authors, and provides evi- 
dence in support of disc regulation, or a similar mechanism, 
operating over the early-PMS stages of the evolution. 

We note that the age of the ONC is still a subject of con- 
troversy in t he literature at the time of writing. S ome recent 
studies (e.g. ISandstrom et al ] |2007l : IJeffriell2007l ) have indi- 
cated that the distance to the ONC may be ~ 3 90 pc, rather 
than t he conventional value of 470 pc used by IHillenbrandl 
l|l997h . Depending on the PMS evolutionary tracks used, 
this can raise the derived age to ~ 2 Myr, which may help 
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Figure 13. Rotation period as a function of mass, using the 
model presented in SI5.2.2I to evolve the ONC distribution (top 
panel) forward in time from 1 Myr to 5 Myr (middle panel), 
and the observed NGC 2362 distribution for comparison (bottom 
panel). Note that these models assume no additional sources of 
angular momentum loss on the early-PMS. 



to explain these apparent discrepancies. We have explored 
this further in Figure IT4l which shows the result of assum- 
ing a 2 Myr age for the ONC, and evolving the rotation 
period distribution forward in time. Since the models then 
predict less spin-up, it is not surprising that the slow ro- 
tators at solar mass are better-fit, but still with evidence 
for (more modest) angular momentum loss being required 
to fully-reproduce the NGC 2362 slow rotator population. 
However, at lower masses, the resulting slower rotation for 
the bulk of these objects places the upper limit to the rota- 
tion periods even higher than observed in NGC 2362. Rais- 
ing the age of the ONC does not therefore readily provide 
an explanation for the problems at the low-mass end that 
have already been discussed. 

We have also repeated the analysis using NGC 2264 as 
the initial condition. Figure [15] shows the results, where we 
have assumed a 2 Myr age for this cluster. NGC 2264 has 
much less differentially-variable reddening than the ONC, 
so the low-mass end of the diagram should be more reliably 
reproduced in this cluster. 

The diagram indicates that this is indeed the case, and 
that the sloping relation for the slow rotators at low masses 
in NGC 2362 is well-reproduced by evolving the NGC 2264 
observations forward in time. Again however, the slowest 
rotators at higher masses (M>0.6 Mq) appear to rotate 
too rapidly compared to the NGC 2362 observations. 

Given that these models do not incorporate additional 
sources of angular momentum loss applying to the early- 




Figure 14. As Figure IT3l except assuming a 2 Myr age for the 
ONC. 



. >. N2264 

: .-. •. • • • ""v^cftft 

• • •• • _ • » • ^ . ^ "W* f • _ " 

: • t *l* r . 



X! 

o 



• • • SB'mode 
» • 



• * •••• •:•*'?.»<•..: ». 



N2362 - 



0.8 



0.6 0.4 

M / M 



0.2 



Figure 15. As Figure H"3l except using NGC 2264 as the initial 
condition, assuming an age of 2 Myr. 
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PMS, such as disc regulation, the ability of the models to 
reproduce the form of the sloping relation for the lowest- 
mass stars (M<0.5 M©) indicates that these stars do not 
experience significant angular momentum loss at this stage 
of their evolution, in contrast to the higher-mass stars. 

We also demonstrated in our ear l ier NGC 2516 and 
NGC 2547 publications ijlrwin et alJ l2007bl : llrwin et all 
|2007cT ). that a similar statement applies to the lowest-mass 
stars at later stages in their evolution: even from the ~ 
5 Myr age of NGC 2362 to the ~ 150 Myr age of NGC 2516, 
additional angular momentum loss at M< 0.5 M© above 
that expected due to (saturated) stellar winds does not ap- 
pear to be required. Given that the higher-mass stars do 
appear to need additional angular momentum losses, this 
implies that such losses are mass-dependent, in the sense 
that lower-mass stars experience less PMS angular momen- 
tum loss, or that they lose angular momentum over a shorter 
timescale (-C 5 Myr) than the higher-mass stars. 



5.3 Comparison with mid-IR Spitzer observations 

In order to search for the correlation between rotation pe- 
riod and the presence of discs expected in the disc reg- 
ulation paradigm for early-PMS angular momentum evo- 
lution, we combined the pres ent optical data-set with th e 
Spitzer IR AC photometry of IDahm fc Hillenbrand! (|2007h . 
Note that IDahm fc Hillenbrand! (120071 1 only give mid-IR 
measurements for objects detected in their combined photo- 
metric and s pectroscopic survey of the central ~ 11' x 11' of 
the cluster (iDahml 12005). whereas the Spitzer field-of-view 
is somewhat larger, so this decreases the available sample 
size. A re-analysis of the Spitzer observations to resolve this 
problem will be attempted in a future publication. 

There are 64 sta r s in c ommon with the catalogue of 
IDahm fc Hillenbrand! l|2007h which have rotation period 
measurements. However, in order to obtain a meaningful 
result, we must restrict the mass range over which the sam- 
ples are compared: not doing so would introduce a bias at 
the faint end toward those objects with excesses, since they 
are more r eadily detected in the mid-I R data. Examining 
Figure 1 of IDahm fc Hillenbrand! (|2007h indicates that the 
I RAC detection l imit i s at [4.5] ~ 15, and using the models 
of iBaraffe et all l|l998h for the M-band (which is a reason- 
able approximation to the Spitzer 4.5 /im band), this corre- 
sponds to a conservative limit of M > 0.4 M©, which has 
been adopted henceforth. 

We have considered two commonly-used indicators of 
the presence of a circumstellar disc: the slope of the IRAC 
spectral ene rgy distribution (SEP ), a = dlog(AF\) / dlog A, 
as used by IDahm fc Hillenbrand! (l2007h. an d the [ 3.6] - 
[8.0] colour, as used by ICieza fc Baliber] (|2006 ) and 
Cieza fc Baliberj (|2007f ). Figure \W\ shows a plot of these 
quantities as a function of rotation period. 

Thresholds were define d in these quantities, following 
IDahm fc Hillenbrand! (|2007f ) for a, with stellar photospheres 
having a < —2.56, "weak discs" —2.56 < a < —1.80, 
and "primordial" or opticall y thick discs a > —1. 80. In 
[3.6] - [8.0] colour, we follow ICieza fc Baliber! (|2007f ), and 
use a threshold of [3.6] — [8.0] > 0.7 to identify objects with 
mid-IR excesses indicative of the presence of discs. Figure 
1161 clearly indicates that the "primordial" discs as defined 
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Figure 16. Two disc indicators plotted as a function of rotation 
period. Symbol sizes indicate the object masses, from 0.4 M© 
(smallest symbols) to 1.2 M© (largest symbols). Top: SED slope 
a as a function of rotation period. Bottom: [3.6] — [8.0] colour 
as a function of rotation period. The dashed lines indicate the 
approximate boundaries in these quantities between regions of the 
diagram occupied by objects with and without discs, as discussed 
in the text. 



by SED slope, and the excess objects in [3.6] — [8.0], are 
essentially the same. 

Both panels of Figure [161 indicate that the 3 — 4 objects 
showing strong disc signatures all have long periods, whereas 
the remainder of the objects show a range of periods, as 
would be expected in the disc regulation paradigm. This is 
also seen in distributions of rotation periods in these bins, 
which are shown in Figure [171 

Applying a two-sided K-S test to the distributions of 
objects with and without IR excesses gave probabilities of 
0.12 (using a > —1.80 to select the objects with significant 
excesses, corresponding to the "primordial disc" objects in 
Figure [T5J , and . 1 5 (using [3 . 6] - [8 . 0] > . 7) that the distri- 
butions were drawn from the same parent population. Since 
the K-S test is not strictly appropriate for such a small sam- 
ple size, we also considered a simple probabilistic argument: 
of the 55 objects without excesses, 30 have periods > 4 days. 
The probability that all three objects with excesses were 
drawn from this same parent population is therefore 0.16. 
These results indicate that the differences are not statisti- 
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Period (days) 

Figure 17. Histograms of rotation period for M > 0.4 Mq, show- 
ing: objects with IR excesses ("primordial discs", or a > —1.80: 
solid line), objects without large IR excesses ("weak discs" and 
"photospheres", or a < —1.80: dashed line), and the entire sample 
for comparison (dotted line). Poisson error bars are shown. 



cally significant. This is most likely due to the small sample 
size. 

In spite of this, the distributions do indeed hint at the 
expected result for the disc regulation paradigm, that the 
objects with discs shoul d be slow rotato r s, as seen e.g. by 
iRebull et all (120061 ) and lCieza fc Baliberl (120071 ). Obtaining 
a larger sample of objects with rotation period and mid- 
IR measurements should improve the significance of these 
results, and allow more definitive conclusions to be made. 



6 CONCLUSIONS 

We have reported on results of an i-band photometric survey 
of NGC 2362, covering ~ 0.4 sq.deg of the cluster. Selection 
of candidate members in a V versus V — I colour- magnitude 
diagram using an empirical fit to the cluster sequence found 
1800 candidate members, over a V magnitude range of 16 < 
V < 26 (covering masses from 1.2 Mq down to below the 
brown dwarf limit). The likely field contamination level was 
estimated using a simulated cat alogue of field obj ects from 
the Besancon Galactic models IjRobin et ajj 120031 ). finding 
an overall contamination level of ~ 65 per cent, implying 
that there are ~ 630 real cluster members over this mass 
range in our field-of-view. 

We derived light curves for ~ 85 000 objects in the NGC 
2362 field, achieving a precision of < 1 per cent per data 
point over 15 < i < 19. The light curves of our candidate 
cluster members were searched for periodic variability corre- 
sponding to stellar rotation, giving 271 detections over the 
mass range 0.1 < M/M < 1.2. 

The rotation period distribution as a function of mass 
was found to show a clear mass-dependent morphology, as 
seen in our earlier NGC 2547, NGC 2516 and M34 stud- 
ies, and qualitatively similar to that in NGC 2264, which is 
thought to have a similar age. 

Using r otational evolution models from our earlier NGC 
2516 work (jlrwin et al 1 l2007bl ). we were able to show that 
distribution of rotation rates for the rapid rotators in NGC 
2362 could be reproduced from the distributions in the ONC 



and NGC 2264 by evolving them forward in time to ~ 5 Myr. 
However, for the slow rotators at masses > 0.6 Mq, the pre- 
dicted rotation rates from the models were too rapid, with 
the observations indicating that the rotation periods of these 
stars are in fact approximately constant over the ~ 1 — 5 Myr 
age range. This confirms the need for additional angular mo- 
mentum losses on this part of the PMS, e.g. as would be 
expected in the disc regulation paradigm. 

Finally, we hav e used the Spitzer observations of 
iDahm fc Hillenbrand! l|2007h to test the disc regulation hy- 
pothesis, by searching for a correlation between rotation rate 
and mid-IR excesses indicative of circumstellar discs. No sta- 
tistically significant correlation was found with the present 
small sample sizes, but the 3 — 4 objects with significant disc 
signatures in the sample are indeed all slow rotators. 
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